The manifestation of quantum resonances and antiresonances 
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We investigate the effect of temperature on resonant and antiresonant dynamics in a dilute atomic gas kicked 
periodically by a standing wave laser field. Our numerical calculations are based on a Monte Carlo method for 
an incoherent mixture of non-interacting plane waves, and show that the atomic dynamics are highly sensitive to 
the initial momentum width of the gas. We explain this sensitivity by examining the time evolution of individual 
atomic centre of mass momentum eigenstates with varying quasimomentum, and we have determined analytic 
expressions for the evolution of the second-order momentum moment to illustrate the range of behaviours. 
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I. INTRODUCTION 

A dilute atomic gas subjected to a pulsed optical standing 
wave is, under certain conditions, an accurate realization of a 
periodically kicked quantum rotor |[T]|2]I2]. The kicked rotor 
is a physical realization of the standard map, which is a well 
established model system in the study of classically chaotic 
Hamiltonian dynamics |4, 5|. The analogous quantum kicked 
rotor provides the opportunity to probe quantum dynamics in 
regimes where the underlying classical system is chaotic |5 

El. 

The atom-optical realization of the kicked rotor, where 
laser-cooled atoms are peiiodically driven by a far-detuned 
laser, is a very attractive system to study because of the pre- 
cise control available in experiments. The amplitude, phase, 
and length of the kicks can be varied over a wide param- 
eter regime, and the initial state of the gas can be care- 
fully chosen using the sophisticated cooling and confine- 
ment methods now available |7|. A large number of such 
laser-driven cold atom experiments have been carried out 

iii]i2i[iii3[ioiiiiiii3[i3iii[i5][i3[i7iii8i[iai2ni 

|2l]|22|23l|24l|25l|26l|22l|2a|23[30l|H 

these have led to the observation of quantum resonant dynam- 
ics IS Em [IT], dynamical localisation H] H [3] [m [HI [HI 
[l4][l5][T6l[T2l[l8l, higher-order resonances lfT9l . and quantum 
accelerator modes ElllQlEIlESIISEHia. In the ultra-cold 
regime, a Bose-Einstein condensed cloud has made it possible 
to observe quantum antiresonances [27 1 and provide detailed 
information about the momentum distribution of a <5-kicked 

gas EHlEa. 

Such atom-optical experimental configurations are clearly 
not exact realizations of the ideal quantum (5-kicked rotor. 
For example, absolutely instantaneous kicks are impossible to 
achieve, and interatomic collisions may also come into play. 
A more fundamental difference is the fact that a laser-driven 
atom models a 5-kicked free particle, rather than a rotor The 
angular momentum spectrum of a quantum rotor is discrete, 
whereas the momentum spectrum of a free particle is contin- 
uous. As first noted by Bharucha et al. this can sub- 
stantially change the manifestation of quantum resonant and 
antiresonant dynamics. 

In this paper we investigate the quantum-resonant and an- 



tiresonant behaviour of clouds of 5-kicked cold atoms | 
El [381 |39l SSI il] Hal |43J as the width of the initial centre 
of mass atomic momentum distribution is varied. In partic- 
ular, we investigate the transition from the zero-temperature 
limit, where the expected behaviour of a (J-kicked rotor is ex- 
actly reproduced, to a qualitatively different behaviour typical 
of atomic samples initially laser-cooled to microKelvin tem- 
peratures. 

The paper is organised as follows. In Sec.|Il]we outline our 
theoretical model for the atom-optical kicked rotor. We begin 



by deriving our model Hamiltonian for the system in Sec. II A 



and give an analytic description of the temporal dynamics in 



Sec. IIB In Sec. II C we outline our numerical implementa- 
tion of the theory, and reproduce some known results to illus- 



trate the method. In Sec. Ill we consider the dynamics of a 



peiiodically kicked plane wave and present a detailed survey 
of the effect of the initial momentum on the system behaviour. 
We derive analytic expressions for the positions and widths 
of the resonances and antiresonances observed, and identify 
a number of higher-order antiresonance features. In Sec. IV 
we model periodically kicked atom clouds with a finite initial 
momentum width. We consider the temperature dependence 
of the system dynamics, and explain the results based on our 
understanding of resonances and antiresonances in the plane- 
wave case. 



II. BACKGROUND 
A. Model Hamiltonian 

1. Two-level atom in a laser field 

We consider a cloud of trapped and laser-cooled thermal 
alkali atoms. The cloud is released from all external fields, 
and then addressed by two far-detuned laser beams that form 
a pulsed sinusoidal potential aligned perpendicular to grav- 
ity |66|. We neglect atomic collisions in the gas, allowing for 
a theoretical description using a single-particle Hamiltonian. 
This is reasonable for a dilute thermal gas. When consider- 
ing experiments carried out on an initially Bose-Einstein con- 
densed sample l27l [28l [29l iH l45|, a noninteracting model 
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serves as a useful starting point. Furthermore, near a Fesh- 
bach resonance, magnetic fields can be used to tune the scat- 
tering length, in principle to zero (e.g., ||451|47]|48]|49|). The 
cloud dynamics are described using a one-dimensional treat- 
ment along the axis of the standing wave. In our noninter- 
acting model, this approach applies without loss of generality 
because the equation of motion for the three spatial degrees of 
freedom is separable. 

Our model Hamiltonian H describes a single two-level 
atom of mass M, with internal ground state \g) and excited 
state \e). The internal atomic levels are coupled by two equal- 
frequency laser travelling waves with a controllable phase dif- 
ference ||50l . As we always consider the lasers to be far- 
detuned (see Sec. |II A2| l, we neglect spontaneous emission 
from the outset. Hence, 



HJ-ncoQ{\e){e\-\g){g\) ■ ^ 



+ -nni \e 
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(1) 



e}{g\+U.c. 



where x is the one-dimensional centre-of-mass atomic posi- 
tion operator with p its conjugate momentum, t is the time, 
Hojo is the energy difference between the internal atomic lev- 
els, Qi and Q.2 are the laser Rabi frequencies with <pi and 02 
the respective laser phases, ojl is the laser frequency, and is 
the laser wave- vector magnitude along the x direction. 

Assuming that the two laser fields are of equal intensity and 
phase (i.e., Q.i - Q.2 - Q./2 and 0i = 02 = </>), 

%2 



H ^^naJo(\e}{e\ - \g){g\) + ^ 



+ -nncosikix) [e^'f'^^'^'^Vx^l + H.c] . 



(2) 



To eliminate the time-dependence in Eq. Q, and shift the en- 
ergy zero to the atomic ground state, we carry out a unitary 
transformation defined by Ui = exp{i[ojL\e}{e\ - a>(){\e}{e\ + 
\g}{g\)/2]t). Hamiltonian Q then transforms to 

H' = nA\eXe\ + + -ftQcos(M) [e'%){g\ + H.c] , (3) 
where we have defined the detuning A = wo - wl. 



2. Adiabatic elimination of the excited state 

In the far-detuned limit (O/A <k 1), and for the case that 
all the atomic population is initially in the internal ground 
state \g), the excited state population is always negligible. In 
this regime, the excited state \e) can be adiabatically elim- 
inated ifsTTl . If we also invoke the unitary transformation 
U2 - exp(-/Q^|g)(g|f/8A), the Hamiltonian simplifies to 



H" 



■ co&iKx), 



(4) 



2M 8A 

where K - 2ki. In Eq. (|4| we consider only the centre-of- 
mass motion of the atoms, as they are always assumed to be 
in the internal ground state. 



3. Time dependence of the laser intensity 

We consider the situation where the laser standing wave 
amplitude, turned on instantaneously at f = 0, is switched 
on and off periodically. Provided that the spectral width 
of each pulse is sufficiently narrow, i.e., the pulse duration 
tp » 1/A, the time-dependent system can be described by 
Hamiltonian Q with an additional time-dependent function 
fit) included in the standing-wave term. The Hamiltonian is 
then 



»2 no} 

Hs= cos(Kx)f(t). 

2M 8A 'J"^' 



(5) 



In the Raman-Nath regime, where the atomic centre-of-mass 
motion is negligible over the pulse duration, /(f) can be ap- 
proximated by a train of (5-kicks |[ll[T3] ll67l . i.e.. 



f{t)^tpYj^{t-nT). 



(6) 



«=o 



In Eq. (|6]l, T is the time between successive pulses, and we 
scale by the pulse duration tp such that /(f) has unit amplitude 



in the limit f „ 



0. 



Finally, our model Hamiltonian for a thermal atomic gas 
subjected to an far-detuned pulsed standing-wave potential 
is L68J 

where we have defined an effective potential depth 0^ = 
Q^fp/8A. Equation ^ is the Hamiltonian for a quantum 6- 
kicked rotor |5 1 (or, more correctly, a (5-kicked particle, as the 
momentum spectrum is continuous rather than discrete). 



4. Spatial periodicity of the Hamiltonian 

Hamiltonian (jTji is spatially periodic. This implies conser- 
vation of the quasimomentum, i.e., Bloch's theorem applies 
1^521. To explain this, it is convenient to expand the position 
and momentum operators into discrete and continuous com- 
ponents: 



X = K-\2ni+e), 
p = HKiic + P). 



(8) 



The eigenvalues of I and k are integers, the eigenvalues of 
9 e [-n, n), and the quasimomentum operator yS is defined 
to take the eigenvalues of the continuous momentum compo- 
nent such that y6 6 [-1/2, 1/2). Substituting the reformulated 
operators ([8]l into Hamiltonian (jTji, yields 



[hKik + I3)f 
2M 



- n^d cos(e) ^ 6{t - nT). 



(9) 



«=() 



The operator /3 commutes with both 6 and k |53|, and hence 
with Hg. Thus, in the system of a (J-kicked particle the laser 
field induces coupling only between momentum eigenstates 
differing in momenta by integer multiples of HK [69] . 
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B. Time evolution in tlie presence of quantum resonances and 
antiresonances 

1. The Floquet operator 

The Floquet operator is the unitary time-evolution operator 
for one temporal period T of the system. The convention we 
have adopted defines the Floquet operator from just before 
one kick to just before the next. The Floquet operator for the 
system governed by Hamiltonian (|9]l is then 



F — exp 



ihK^T. 
2M 



cos 9). 



(10) 



We consider F applied to a state with a particular value of the 
quasimomentum /3, i.e.. 



k=-oo 



(11) 



Note that the dimensionless momentum eigenket \k + (3) satis- 
fies the eigenvalue equation 



p\k +/3}^ nK(k + /3)\k + P), 
and the orthogonality condition 

{k' +l}'\k+p)^6kk'6{Ji-li'). 



(12) 



(13) 



As the q uasimomentum is a conserved quantity (see 
Sec. nA4i, the operator p in the Floquet operator (lOi can 
be replaced by its eigenvalue /? when acting on the quasimo- 



mentum eigenstate (111, i.e., 

Fi<A(y6)) = Fmm), 

where 

, ihK^T 

F(J5) = exp 



(14) 



2M 



[/t -H /?]^ exp(/0rf cos 0). (15) 



2. Quantum resonances and antiresonances 

The time evolution described by the Floquet operator ([TSjl 
can vary significantly depending on the system parameters. 
We consider the strong driving regime, where the classical 
stochasticity parameter k > 1.5 1 6 1. In terms of parame- 
ters relevant to the quantum system described in this paper, 
K = K^Tticpd/M. In the case of strong driving dynamical lo- 
calization is typically observed ia]|2l|3][n][I2[l3l[ll[l5l 
[16] [nl [18] |54] |55] |56] |52l, whereby the atomic momentum 
spread increases linearly with time at the classical "diffusion" 
rate [4], i.e., 



Kkt) Y 



n + </3 >o, 



(16) 



up until a quantum break time jTJ [121 155] . While dynamical 
localization is the more common scenario, the Floquet oper- 
ator simplifies considerably for particular choices of the time 



T between successive 5-kicks. Distinctive resonances and an- 
tiresonances then occur in the 5-kicked particle system (e.g., 
|[I][8][l9l|27]|29l). 

To illustrate resonance and antiresonance in an atom-optical 
context we consider their simplest manifestation, which oc- 
curs in the fi - Q subspace; this means assuming the initial 
system state |i^(/? = 0)) [see Eq. ([TT])]. Due to quasimomen- 
tum conservation, fi must always remain = 0. The accessible 
momentum spectrum is therefore discrete, consisting of inte- 
ger multiples of TiK only. This can be mapped to the spectrum 
of an angular momentum, which consists of integer multiples 
of h.. The results presented here for the fi - Q subspace are 
therefore equivalent to those gained from theoretical quantum 
studies of the "genuine" rotor |3] |58]| . 

We consider the (5-kicked particle for the case where the 
pulse periodicity T is given by 



{2nM 



(17) 



where ^ is a positive integer. The Floquet operator ( [TSj l, with 
yS = 0, is then F{Q) - exp(-i/'7rP) exp(/0^ cos 0). The free 
evolution operator exp(-;^7r^^) is equivalent to txp{-i(nk) - 
expiifnk), because the eigenvalues of k are integer, and parity 
is conserved when an integer is squared. The time evolution 
due to the Floquet operator F(Q) - expi+Knk) exp(/0j cos 6) 
is dramatically different depending on whether i is even or 
odd. 

For the case that ( is even, the free evolution operator col- 
lapses to the identity, i.e., F(Q) = exp(i(pj cos 9), and 



F"(Q) = e 



irKpd cos ^ 



(18) 



Thus, after a time riT, the final system state is as if a single 
<5-kick had been applied with a strength n times the original 
kick strength (pd- This phenomenon is known as a quantum 
resonance. The shortest pulse periodicity for which resonance 
occurs, i.e., Tj - AnMjTiK^ {£ - 2), is known as the Talbot 
time by analogy with the Talbot length in classical optics 1591 . 

Considering now the general case where I may be odd or 
even, we note that [9, k} - i ||53]| , and therefore obtain 



(19) 



If £ is even, we regain Eq. ( 18 i for n = 2. If ^ is odd, 
F^(Q) = 1, and the system returns to its initial state every 
second pulse. The recurrence of the initial state with period 
2T is the phenomenon known as antiresonance. Higher-order 
antiresonances are also possible, where the initial state recurs 
after a larger number of iterations of the Floquet operator (see 
Sec.jm]) |70|. 

In this paper we consider the 5-kicked particle for pulse 
periodicities T - {Tt/2, where £ is a positive integer [see 
Eq. ( 17 1]. For the /3 - subspace, resonance occurs for even 
values of and antiresonance for odd values of {. However, 



for a general p subspace, the Floquet operator ( 15 i is 



F(J3) = e-'^P V 



OT/J- C -m(k{ 1 +213) '4 J cos 6 



(20) 
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The time evolution of the ^-kicked particle, as governed by 
the Floquet operator pO]), can be evaluated analytically. In 



Sees. IIB3 and II B 4| we derive expressions for the system 



evolution, which encompass results presented previously by 
Bharucha et al. Ill 111 for even {. 



3. Momentum eigenstate evolution 



In this section we evaluate the time-evolution of a <5-kicked 



particle, as described by the Floquet operator (20 1. We con- 
sider a system initially prepared in the momentum eigenstate 



m^o))^\k+i3). 



(21) 



Note that, once the evolution of state ( |2T] i is known, the evo- 
lution of any quasimomentum eigenstate [i.e., Eq. ([TT]l], and 
in fact any general state, can also be determined. A general 



treatment will be presented in Sec. II B 4 



After n iterations of the Floquet operator, state \^{t - 0)) 
becomes 



mt^nT))^F"(l3)\k+p). 



(22) 



Due to quasimomentum conservation (see Sec. |II A 4| i, 
Eq. ( 22 1 can be expanded using momentum eigenstates with 



the initial quasimomentum fi, i.e., 

oo 



(23) 



j=- 



where the probability amplitudes Q;(/6, nT) are given by 

Ckj{p,nT)6(p-p') = {j+p'\F"mk+li)- (24) 

We impose the normalisation condition 2y |Qy(y6, «r)p - 1. 

The complex amplitudes ciij{[5,nT) can be most simply 
evaluated using Eq. ( 24 1 and 



{j+/3'\F"(m+f3}^ jdx{j+/3'\x}{x\F"(m+/^\ (25) 
where 



(26) 



A lengthy deduction (see Appendix [A]) dete rmines that the 
matrix element {x\F"(fi)\k + /3} [see Eq. \A21\ ] is 



{x\F"mk+/S}=- 



ilnpr -innp-t i(k+l3)(Kx-2nT) 



sin(nY) 
sinY 



J^-'V("+l)T 



(27) 



where, for convenience, we have defined the jS-dependent vari- 
able 



T= l;r(l+2;S)^. 



(28) 



The matrix element ( [25] t then becomes 
{j + fS'\F"(J3)\k + p) =e-''"'/*'fe-'2«*if 



sin(nT)'^ 



X I 



sinT / 
27r J 

(29) 



Evaluating the Fourier integral in E q. (|29] l [see Eq. ( 13 1], and 
substituting the result back into Eq. (|24|i, we find that 

ta T\ T [a. sin(nT)\ -,(/-yfc)(«+l)T -Imtk -iimB^e 

Ckj{p,nT) ^ J j-k\4>d ■ ^ M e " " ' e e . 
\ smT / 

(30) 

For the system initially prepared in the momentum eigen- 
state pT) , the ^th-order momentum moment after time nT is 



{p%, = (hKf 2 \cujifi,nT)\\j + pf 



j=- 



Jj-k 4>d 



sin(«T) 
sinY 



(31) 



U+P)"- 



The distinctiveness of resonant and antiresonant behaviour in 
the 5-kicked particle can be concisely characterised by con- 
sidering the time evolution of the second-order momentum 



moment [Eq. (31 1 with q - 2], which is proportional to the 
kinetic energy. 



4. Evolution of an incoherent mixture 

The above treatment can be generaUsed to the case of an 
atom cloud at finite temperature, by considering an incoherent 
mixture of plane waves with an initial momentum distribution 
D(nK(k+f3)) = Dk(/3)/nK. The initial density operator for the 
system is taken to be 



p^nKj dp\p}Dip){p\, 

-1/2 



dp y \k+p)Dk(mk+p\, 

1/2 kt-L 



(32) 



where we have assumed that the density matrix in the mo- 
mentum representation is initially diagonal, i.e., there are no 
coherences in the system (as is the case for a thermal state). 
After time nT, the density operator evolves to 



Xl/2 oo oo 

dpY, J]c,j(p,nT)\j+p} 
1/2 k=-ca [j=-ca 



DkiP) 



£ cl.,{J3,nT){j' + p\ 



(33) 



where cj^jiPjiT) is given by Eq. (30i. Considering the diag- 
onal elements of the density operator (33i, we find that the 
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momentum distribution evolves to 

DO 

D,{J3, t^nT)^Yj lO't*^' nTfDjdS), (34) 



j=- 



and the g'th-order momentum moment is given by 



dp V \cuj{l3,nT)\ Dumi+Pf- (35) 

1/2 ).*t^co 



C. Numerical implementation 

1. Monte-Carlo method 

The numerical results presented in this paper are generated 
using a Monte Carlo approach. This is used in preference to 
a grid-based method, which could produce sampling errors if, 
for example, a disproportionate number of grid points were to 
coincide with strongly resonant or antiresonant values of the 
initial momentum. 

The initial condition consists of N plane waves distributed 



according to the momentum distribution Dii(J3) (see Sec. II B 4 



and Sec. II C 3 1. The time evolution of each momentum eigen- 



state is evaluated by applying the Floquet operator ( 20 1, which 
is represented in a momentum basis |9|. The final momen- 



tum distribution [see Eq. ( 34 1] is constructed by averaging the 
distributions resulting from each of the individual plane-wave 
time evolutions. 

To illustrate our method we consider resonance and antires- 



onance for (i) the yS = subspace (as described in Sec. 11 B 2 1, 
and (ii) a cold (but thermal) atomic gas [9|. Sections II C 2 



and II C 3 provide an introduction to the main content of this 
paper by summarising some relevant well known results (e.g., 
ILiSJ). 



2. Resonance and antiresonance in the /? = subspace 



In Sec. II B 2 we showed how the dynamics of the <5-kicked 
particle depends dramatically on the pulse periodicity T . In 
particular, resonance occurs in the /? = subspace for even £ 
(where T - {TtI2 [see Eq. ([iTjl]), and antiresonance is ob- 
served for odd values of {. Figure[T]shows the atomic momen- 
tum distribution and the second-order momentum moment of 
a 5-kicked particle, for the case where the system is initially in 
the zero-momentum eigenstate with unit amplitude. We have 
chosen (p^ - 0.8;r as an illustrative value typical of recent ex- 
periments ||3]|9l|20l|2T]|22l|2l|Mll25l|26l. 

Figures [ija) and (b) are for the case ( — \. We observe the 
characteristic feature of an antiresonance, i.e., recurrence of 
the initial state with period 2T . In contrast. Figs. [TJe) and (f) 
correspond to the resonant case { - 2. We observe quadratic 
growth in the second-order momentum moment, otherwise re- 
ferred to as ballistic expansion |8|. Figures [TJc) and (d) show 
the system evolution for the case where T is an irrational mul- 
tiple of the Talbot time. We observe dynamical localisation. 
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FIG. 1: Momentum distributions [(a), (c), and (e)] and correspond- 
ing second-order momentum moments [(b), (d), and (f)] of a (5-kicked 
particle. In each case the initial condition is a zero-momentum eigen- 
state. Parameters are N = I, (pj = 0.8n, and (a), (b) T = Tx/2, (c), 
(d) r = (1 + V5)rr/2, and (e), (f) r = Tr. The solid lines in (b), 



16 1 (with (POo = 0), 



(d), and (f) are analytic results from Eqs. 1 54 1, 
and \50\, respectively. 



characterised by classical-like diffusion until a quantum break 

time Em mini Ea. 

As described in Sec. |IIB2 due to quasimomentum conser- 
vation the accessible momentum spectrum for the dynamics 
considered in Fig.[T]can be mapped to the spectrum of an an- 
gular momentum. The results presented in Fig. [T] are there- 
fore equivalent to those that would be obtained by studying 
the quantum (5-kicked dynamics of a "genuine" rotor Il5l l58ll . 



3. Resonance phenomena for a finite temperature cloud 

In the case of a finite temperature cloud, all the quasimo- 
mentum subspaces are initially populated, not just the /3 - 
subspace. Note that the Floquet operator governing the dy- 
namics in each quasimomentum subspace is explicitly depen- 



dent on /3 [see Eq. ( 15 i]. The overall dynamics can therefore 
exhibit some significant qualitative differences when com- 
pared to the case where the dynamics are restricted to the 
j0 = subspace (and, by extension, the dynamics of a "gen- 
uine" rotor). 

We model a cold thermal cloud of atoms by considering a 
Gaussian initial momentum distribution Ditifi) with zero mean 
and standard deviation w: 



1 



■V2^ 



exp 



2w2 



(36) 



This corresponds to a Maxwell-Boltzmann distribution for 
free particles with temperature 7~„, - h^K^w^ /rnks- The 
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second-order momentum moment of the kicked cloud, after 
time riT, is determined by substituting the momentum distri- 



bution ( 36 1 into the general momentum moment evolution ex- 
pression ( 35 1 for ^ = 2, to give 



pl/2 - 
-'-1/2 J,t^o. 



W 



'j-k 



sin(nT) 
' sinT 



(37) 



Figure |2] shows the atomic momentum distribution and the 
second-order momentum moment of a 5-kicked cloud of non- 
interacting particles with initial standard deviation w - 2.5 
(corresponding to approximately 5 ;uK in the case of caesium 
[9 1). In contrast to the zero-momentum eigenstate case (see 
Fig.[T]i, the finite temperature cloud appears to behave identi- 
cally for ( = 1 [Figs, ga) and (b)] and f = 2 [Figs, [^e) and 
(f)] |71J. For £ - 1 or 2, a small fraction of the cloud expands 
ballistically, but much of the distribution remains near zero 
momentum. This clustering near the centre is due to antireso- 
nant effects occuring concurrently with the quantum resonant 
dynamics. The result is limited ballistic expansion, which will 
be discussed in more detail in Sec. |IV B] In Figs.|2|c) and (d), 
where T is an irrational multiple of the Talbot time, we ob- 
serve dynamical localization rather than ballistic expansion. 
The distribution in Fig.[2|c), while truncated at the edges, ap- 
pears flatter and broader compared to Figs.|2|a) and (e). 




10 20 30 

n 



FIG. 2: Momentum distributions [(a), (c), and (e)] and corresponding 
second-order momentum moments [(b), (d), and (f)] of a 5-kicked 
particle. In each case the initial atomic momentum distribution is 
given by Eq. l[36j. Parameters are N = 10000, (fid = 0.8;r, w = 2.5, 
and (a), (b) T = Tt/2, (c), (d) T = (1 + V5)rr/2, and (e), (fjT = Tt. 
The solid lines in (b), (d), and (f) correspond to the linear result ( |16^ . 

Figures [2|b) and (f) show that the second-order momen- 
tum moment grows linearly. This is dramatically diff'erent 
from either the periodic behaviour or the quadratic growth 



observed in Figs. [TJb) and (f), respectively. Linear energy 
growth is characteristic of resonant phenomena in a thermal 
atomic cloud |[8l, and the slope of the line [see Eq. ( 16 1] is 
well approximated by the classical diffusion rate 1 1 , 4|. Close 
inspection of Figs. |2|b) and (f) reveals that the classical ex- 
pression slightly overestimates the numerical data in the an- 
tiresonant case, and underestimates in the resonant case. This 
will be clarified in Sec.lIVI 



III. 



MOMENTUM DEPENDENCE OF RESONANT AND 
ANTIRESONANT FEATURES 

A. Motivation 



In Sees. II B 2 and II C 2 we considered the j6 = subspace. 



and found that resonance and antiresonance are clearly dis- 
tinct. In contrast, at finite temperature (see Sec. II C 3 i, the 



system appears to behave identically irrespective of whether £ 
is even or odd (T = {Tj/l). In the remainder of this paper we 
investigate in detail the transition between the behaviour ob- 
served in Fig.fTlfor a zero-momentum eigenstate, and that ob- 
served in Fig.|2]for a broad thermal distribution. We consider 
only the cases where the pulse periodicity T is a half-integer 
multiple of the Talbot time (i.e., T = £Tt/2). In this section 
we begin by looking at the effect of the quasimomentum /3 on 
resonant and antiresonant dynamics in the system. 



B. Resonance and antiresonance conditions 

We consider a system initially prepared in a momentum 
eigenstate [see Eq. ([2T]i]. A resonance in the system oc- 
curs when the free evolution operator collapses to the iden- 
tity. From the Floquet operator ( 20 1, ignoring the global /3- 



dependent phase, we find that resonance occurs when 



1 



-(1 +2fi)£ = m, 



(38) 



where m is an integer In general, the probability that the ini- 
tial state \k + P) is measured to have momentum tiKij + /?), 
after a time nT, is determined from Eq. (|30]l to be ifTTIl 



\ckj(fi,nTf 



Jj-k\^d 



sin(nY) 
sinY 



(39) 



At resonance, when Eq. ( 38 1 is satisfied, Y - mn [see 



Eq. (28 1], so sin(nY) = and sinY - 0. Taking the limit 
Y — > mn, Eq. ( 39 1 becomes 



\ck0,nTf = [Jj-k(n<Pd)f. 



(40) 



As discussed in Sec. |II B 2| the system responds as if a single 
<5-kick had been applied with strength ncp^. 



From Eq. ( 39 1 we see that antiresonance occurs when 



^(1 + 2yS)^ = m+ ^, 



(41) 



7 



for integer m. In this case T = (m + 1 /2)jt, and the probability 
that the initial state is measured to have momentum ftK(j + P) 
is 



\cuj(J3,nTf = [y,_,(0rfSin(«7r/2))]2 



(42) 



For even values of n, \ckj(J3, nT)\ - 5jk, and the initial state re- 
curs. Hence, every kick-induced diffraction is exactly inverted 
by a kick-induced contraction. 



From Eq. ( 39 1 it is also apparent that higher-order antires- 
onances are possible. An Mh order antiresonance, for which 
the system returns to its initial state every Mh kick, occurs 
when 



2 ^ 



(43) 



where m is a nonzero integer, and m and have no common 
factors. In this case. 



\cki{J3,nT)\' = 



Jj-k\<Pd 



sininmn/N) 
sm{mn/N) 



When n is a multiple of the order A^, \ckj(J3, nT)\^ 
the system returns to its initial state ll72l . 



C. Role of the quasimomentum 



(44) 



Sjk, and 



Figure [3]illustrates the j3 dependence of the system dynam- 
ics. In Fig.^a), for i - 1, an antiresonance occurs at;8 = 0, as 
observed in Figs.[TJa) and (b). Resonances occur at/5 = +1/2. 
In Fig. |3jb), for ^ = 2, a resonance occurs at y6 = 0, as ob- 
served in Figs.[TJe) and (f), and also at y6 = +1/2. Antireso- 
nances are observed at/3 - ±1/4. 



(a) 



(b) 



I 



I — I 

1100 
50 







FIG. 3: {(/3^>>,'/^ of a (5-kicked particle initially in the plane-wave 
state ^>. Parameters are N = \,4'd = O.Sin', and (a) T = TtI2 and (b) 
T = Tj. 



The yS dependence of the system evolution can be investi- 
gated in more detail by looking at the evolution of momentum 
eigenstates \k + /3), where k - Q and /3 e [-1/2, 1/2). It is 
convenient to consider the momentum variance 



{P)l- 



(45) 



The momentum variance (45 1 displays the same signatures of 
resonant and antiresonant dynamics that have been described 
for the second-order momentum moment [see Figs.[T|b) and 
(f), and Sec. |IIB 3j . However, the variance is more useful in 
the respect that its evolution is independent of the discrete mo- 
mentum k [as we show in Eq. (|46]) below]. In Figs.[T]and|2] the 
mean momentum {p)„ is zero (or negligible) and it is not nec- 
essary to distinguish between the variance ( |45] l and the mo- 
mentum moment {p^)„. However, in this section we consider 
the momentum variance ( |45] l, because the mean momentum 
of the eigenstate \k + /3} is in general nonzero. 

The evolution of the momentum variance for a <5-kicked 
particle initially in a momentum eigenstate \k + /i) can be de- 
termined from Eq. (31 1. Changing the summation index to 
/ = j - k, and using Eqs. (B6 1 and (B7 1, the evolution of the 
variance simplifies to 



sin(nT) 
sinT 



(/ +/3f 



(46) 

The variance (46 1 is independent of k so we can explore the 
evolution of all momentum states \p} by considering only k - 
eigenstates where the quasimomentum /3 spans the range yS e 
[-1/2,1/2). 



The values of /3 for which resonances and antiresonances 
occur can be determined by rearranging the resonance and 
antiresonance conditions, (38i and (41 1, respectively. A res- 



onance occurs in a particular y6 subspace if /3 = /^f^, where 
- m/{ - 1 /2 and m is an integer. The spacing between res- 
onances is then ySjJ^j -/Sf, - l/i. An antiresonance occurs in a 
particular/? subspace if /? = /3^„ where = (2m+ 1)/2£- 1 /2. 
Antiresonances occur exactly halfway between resonances, 
and are also separated by /3^^^ - f^m - 1/^- 



D. Resonance and antiresonance width 

1. Reconstruction loci of higher-order antiresonances 

The resonances and antiresonances in Fig. |3] have a mo- 
mentum width associated with them, which depends on ( and 
the kick number n. As a precursor to deriving expressions for 
these widths, we zoom in on resonant and antiresonant fea- 
tures for £ - \ (see Fig. |4| and 1-2 (see Fig. |5]l. In close 
proximity to the main resonant and antiresonant features, we 
observe hyperbolic curves of zero momentum variance. We 
call these reconstruction loci, as they each correspond to a 
discontinuous series of points in time where higher-order an- 
tiresonances cause the initial condition to be periodically re- 
constructed. 
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It is now possible to define a resonance width as follows. 
The reconstruction loci closest to a resonance are described 
by ;0 = +1/N{. For \m'\ - 1, N and are equivalent, and 
the reconstruction loci can also be written as ^0 = + 1 /NL The 
hyperbolic curves lying halfway between the resonance at yS,^ 
and the reconstruction loci /3 - +1/N{ we call transition loci, 
as they are the sequence of points at which the wave-function 
begins the process of reconstruction. The transition loci are 
described by y6 = +l/2N{. We quantify the resonance width 
by the separation in /3 between the two transition loci adjacent 
to the resonance, i.e.. 



N{ 



(48) 



FIG. 4: {{p^))U^ of a (5-kicked particle initially in the plane-wave 
state Ip). Parameters are N = I, ipd = 0-Sn, and T = Tt/2. 




0.25 

p[nK] 



FIG. 5: ((p^))]J^ of a (5-kicked particle initially in the plane-wave 
state Ip). Parameters are N = I, (f>j = O.Sn, and T = Tj. 



2. Resonance width 

Here we determine the locations of the higher-order an- 
tiresonances, at which the system returns to its initial state 
i{{p^))n^ = 0). We begin by rewriting Eq. as 



1 



(1 +2P){ ^m + m IN, 



(47) 



where m and m' are integers, and < \m'\ < r^/21. The 
integer .^V can be any kick number for which the initial state 
reconstructs as a result of an Mh-order antiresonance, i.e., 
is an integer multiple of the antiresonance order A^. A direct 
correspondence can now be made between the integer m in 
Eq. (47 1 and the integer m in y6^, [or, equivalently, in the res- 
onance condition (38 1]. In close proximity to a resonance, 
we evaluate the difference /3 - /3 - /3^, for which /3 satisfies 
the higher-order antiresonance condition (43 i. We find that 
the terms involving m cancel, and the equations /3 = m' /Ni, 
where < |m'| < r^/21, describe the reconstruction loci in /3 
and n. 



3. Antiresonance width 

A measure of the second-order antiresonance width can 
be determined in a similar way to the resonance case, by 
considering the closest antiresonances (of any order) near an 
N - 2 antiresonance. In general, the higher-order antireso- 
nances in close proximity to an = 2 antiresonance occur for 
P - P - fi'^„, where fi satisfies the higher-order antiresonance 



condition (43 1. It is convenient to again use Eq (47 1, but now 
m' is constrained to the values < m' < N. Making a di- 



rect correspondence between the integer m in Eq. (47 1 and the 
integer m in y6^, we find that yS - (2m' - N)/2N{. Defining 
s = 2m' - N yields p = s/2Ne, where s = (-|A^ - 2|, -|A^ - 
2| -h2,-|A^-2| -h4,...,A?-2). 

The reconstruction loci are described by j6 = s' /2N{, within 
the limits < < ./V, where .^V and s' must be either both 
even, or both odd. Looking closely at the reconstruction loci 
adjacent to the antiresonances in Figs. [4] and [5] it is possi- 
ble to confirm that adjacent reconstruction loci correspond to 
.^V values of opposite parity |73|. In particular, the N - 2 an- 
tiresonance has even parity and the closest reconstruction loci, 
described by ^0 = +\I2NC (taking i' = +1 and A^ = N), have 
odd parity. 

Transition loci can be defined which lie halfway between 
the N - 2 antiresonance at fi^ and the closest reconstruction 
loci described by y6 = +l/2N{. In the antiresonance case the 
transition loci /3 - +l/4N{ are the sequence of points for 
which the even parity of the N - 2 antiresonance changes 
over to the odd parity of the closest higher-order antireso- 
nances. The width of a second-order antiresonance can be 
characterised by the separation in y6 between the transition loci 
P = +\IAN{, i.e.. 



= We 



(49) 



The antiresonance width ( |49| ) is a factor of two narrower than 
the resonance width 6Pr, as is evident from Figs.|4]and|5] 
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IV. EFFECT OF THE INITIAL ATOM CLOUD 
TEMPERATURE 

A. Overview 

In Sec. Hn] we found that the evolution of the momentum 
variance for momentum eigenstates \k + (3) depends dramat- 
ically on the quasimomentum fi, but is independent of the 
discrete momentum k. This means that the evolution of a 6- 
kicked particle, with a Gaussian initial momentum distribu- 
tion [see Eq. ([36]l], depends sensitively on the initial popula- 
tion of /? subspaces B2l . In particular, if the initial momentum 
distribution is broad compared to the momentum separation 
IjE between consecutive resonances (or antiresonances) [see 
Sec. |IIIC| , the whole range of individual momentum eigen- 
state evolutions is covered (from completely resonant to com- 
pletely antiresonant). In this case an averaged evolution re- 
sults [8, 11, 16 1. 

In this section we investigate in detail the effect of the ini- 
tial momentum distribution width w on the evolution of a 5- 
kicked particle, or ensemble of particles. We give analytic 
expressions for the second-order momentum moment evolu- 
tion for both a zero-momentum eigenstate (see Fig. [T]l and for 
a finite temperature cloud (see Fig.|2]). We then use our knowl- 
edge of the p dependence of momentum-eigenstate evolutions 



(described in Sec. Ill i, to understand the changeover between 
the zero-temperature regime where resonant and antiresonant 
features are distinct, and the finite temperature case where the 
system energy is observed to grow linearly [see Figs.|2|b) and 
(f)]. 

Note that for the remainder of this paper we consider atomic 
momentum distributions symmetric around zero. Thus, it is 
equivalent to consider either the second-order momentum mo- 
ment or the momentum variance [see Eq. ([45]l]. We consider 
the former for notational convenience. 



B. ElTect of temperature on resonant dynamics 

For a zero-momentum eigenstate, resonance occurs for 
even I {T - {TtI2). At a resonance, the second-order mo- 
mentum moment evolves quadratically as 



1 



2 2 



(50) 



where we have taken Eq. (B12 1 with /3 - and f even. The 



quadratic dependence of the second-order momentum mo- 
ment with the kick number n is illustrated in Fig [Tff). In 
contrast, the energy of a finite temperature cloud grows lin- 
early at a resonance, as shown in Fig. |2|f). In the latter case, 
the growth rate is well approximated by the classical diff'usion 
result [given by Eq. ([T6]l]. 

Figure|6](symbols) shows the evolution of the second-order 
momentum moment for three initial Gaussian momentum dis- 
tributions with different values of the standard deviation w. 
We observe that the momentum moment initially follows the 
quadratic growth of the zero-momentum eigenstate case (up- 
per dotted line in Fig.|6]l but, for a finite temperature cloud, the 



ballistic expansion becomes limited at a particular kick num- 
ber, denoted ««. Above the rate of energy growth begins 
to revert to the classical diffiision result (16 1 (indicated by the 
lower dotted line in Fig.|6]l. For large w, the energy growth is 
bounded below by the classical result (16 1 [see Fig.|2|f)]. 




FIG. 6: Second-order momentum moment of a ^-kicked atom cloud. 
The initial state is a Gaussian momentum distribution with standard 
deviation (x) w = 1/128, (o) w = 1/32, and (+) w = 1/8. Param- 
eters are N = 10000, if>j = O.Sn, and T = Tj. The solid lines are 
analytic results for square distributions [see Eq. \52\], chosen with 
an initial standard deviation in agreement with the corresponding 
Gaussian case. The upper dotted line corresponds to the quadratic 
growth \5Q) , and the lower dotted line corresponds to the classical 
linear growth \16\ , taking (p^>o = 0. The vertical dashed lines indi- 
cate n = njt [see Eq. ^5l\]. 

The behaviour observed in Fig. |6] can be understood by 
considering the range of /? subspaces involved in the system 
evolution (see Sec. [Illji. The resonance width was found in 
Sec[mD]to be S^r = 1/Nt [see Eq. (HHll]. Provided that the 



initial momentum width of the cloud is well within the reso- 
nance width, resonant dynamics dominate the system evolu- 
tion. Conversely, if the initial momentum width is larger than 
the resonance width, higher-order antiresonances play a sig- 
nificant role and the ballistic expansion is limited. 

Quantitatively, we define a kick number to indicate the 
number of kicks at which the ballistic expansion becomes lim- 
ited. We define by requiring that two standard deviations 
of the initial Gaussian momentum distribution lie within the 
resonance width 6(5r, i.e.. 



1 

4we' 



(51) 



Equation (jSTJi ensures that 95.4 % of the atom cloud expands 
ballistically up until the kick number hr (indicated by the ver- 
tical dashed lines in Fig.|6]l. 

To derive an analytic expression for the second-order mo- 
mentum moment evolution of a finite temperature cloud, we 
consider a <5-kicked particle with a square initial momentum 
distribution (see Appendix |B 2\ . We find that the second-order 
momentum moment evolves as 



■ sin(2m7r^e) 



m=l 



(52) 
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where e is the half-width of the initial momentum distribution, 
and m is integer [see Eq. ( |B18| l with even €\. If the bound- 
aries of the square distribution lie exactly on resonant or an- 
tiresonant values of the quasimomentum (i.e., e = m' I2(, for 
positive integer m' within the limits < m' < {), all the /3 
subspaces are populated equally (see Sec. |I1IC| and Fig. [3]l. 
Equation ( 52 1 then simplifies to 



(53) 



Noting that {p)l = we recover the linear growth in 

the momentum moment evolution, as predicted by the classi- 
cal result ( [T6] l. 



The moment evolution (52 1 of the square distribution is in 



good agreement with the moment evolution of a Gaussian mo- 
mentum distribution with the same initial standard deviation, 
as shown in Fig.|6] The moment evolution of a square distribu- 
tion slightly underestimates the moment evolution for the cor- 
responding Gaussian distribution, because the Gaussian has a 
relatively larger proportion of the initial population in close 
proximity to the resonance (p = 0). 



C. Effect of temperature on antiresonant dynamics 




FIG. 7: Second-order momentum moment of a (J-kicked atom cloud. 
The initial state is a Gaussian momentum distribution with standard 
deviation (a) (□) w = 0, (x) w = 1/128, and (o) w = 1/32, and (b) 
(+) w = 1/8, (*) w = 1/4, (>) w = 1/2, and (o) w = 2.5. Parameters 
are N = 10000, cpj = OM, and T = Tt/2. In (a) the solid lines ai'e 
given by Eq. l |56[ ) and the vertical dashed lines indicate n = [see 
Eq. ^55\]. The solid lines in (b) are given by Eq. \16\ . 



For a zero-momentum eigenstate, antiresonance occurs for 
odd integers {. At the antiresonance, |co)(0)|^ = Jj(<l>d) for odd 
n, and |c(y(0)p = djo for even n [see Eq. (42l]. The second- 
order momentum moment evolves as 



1 



(54) 



n'K'-cl,',[l-(-ir], 



where we have taken Eq. ( |B12| i with y6 = and i odd. The os- 
cillation of the second-order momentum moment with period 
2T is illustrated in Fig[TJb). 

The effect of the initial atomic momentum width on the 
second-order momentum moment evolution at an antireso- 
nance is shown in Fig. |7] The evolution can be separated into 
two regimes: (i) when w <s; 1/2^ periodic oscillations are ob- 
served [see Fig.[7|a)], and (ii) for w > l/2{ the momentum 
moment grows linearly [see Fig. [7|b)]. For large w, the en- 
ergy growth is bounded above by the classical result ( [T6] l [see 
Fig.||b)]. 

As in the resonance case, we can define a kick number 
below which the majority of plane-wave components in the 
distribution antiresonate with period 2T. Requiring that two 
standard deviations of the initial Gaussian atomic momen- 
tum distribution lie within the antiresonance width 6/3a [see 



Eq. (49 1], we find that 



HA = 



(55) 



In Fig.[7|a), the kick numbers are indicated by the vertical 
dashed lines. For kick numbers less than the initial state 
recurrence with period 2T is clearly observed. 



In the low temperature regime [where w l/2{, as in 
Fig. |7ja)], we observe oscillations with period 2T even af- 
ter n - Ha- However, the oscillations decay to a constant 
value. The momentum moment evolution in the low tempera- 
ture regime is well approximated by 



-n K 



tp-KW, (56) 



as shown in Figj7[a). The derivation of Eq. (56 1 is presented 
in Appendix |B 3| 

The observed decay in the periodic oscillations occurs 
when the initial momentum width of the cloud is compara- 
ble to, or larger than, the antiresonance width 5[iA- Provided 
that w <K 1 121, higher-order antiresonances begin to play a 
role in the system dynamics, but resonances do not. In this 
regime the evolution of the different [} subspaces can be de- 
scribed as an oscillation of the same amplitude as the j5 - Q 
subspace, but with a different phase. This means that the dif- 
ferent p subspaces gradually become out of phase with each 
other, leading to a dephasing in the oscillation of the second- 
order momentum moment. The decay of the oscillations has 
been observed experimentally by Duffy et al. |27 1. 

For the case where w > 1 12{ [see Fig.|7]^b)], certain initially 
populated fi subspaces fulfill the resonance condition (38 1, so 



the second-order momentum moment is observed to grow lin- 
early. For sufficiently broad initial momentum distributions, 
the slope of the linear moment growth tends to the classical 



diffusion rate [see Eq. (16 1] 
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V. DISCUSSION AND CONCLUSIONS 

In this paper we have investigated theoretically the quantum 
resonant and antiresonant dynamics of 5-kicked cold atoms. 
In particular we have investigated how the behaviour changes 
as one considers progressively broader initial momentum dis- 
tributions. That is, as one considers increasing the initial 
temperature of the atomic cloud from the ideal case of zero 
Kelvin. We have identified transition times, as a function of 
the initial momentum width, after which the quantum reso- 
nant and antiresonant dynamics clearly deviate from the ideal 
zero temperature case. This can be understood by examin- 
ing the time-evolution of individual momentum eigenstates, 
which varies greatly with the quasimomentum. Very cold tem- 
peratures are typically achieved with an atomic Bose-Einstein 
condensate, which introduces a mean field. For the results 
discussed by Duffy et al. [271, mean field effects are not very 
significant compared to the eff'ect of the sample's initial mo- 
mentum width. The analysis presented in this paper can then 
still be applied. It should be noted, however, that with in- 
creased strength of the interatomic interactions, this may no 
longer hold [601 ISD. 

Typical temperatures achieved with an initially magneto- 
optically trapped sample, such as used in the Oxford exper- 
iment [3, 9, 20, 21. 221 mm |25] Ea (temperature ~ 5 
//K, corresponding to a Gaussian momentum distribution with 
w ~ 2.5hK) observe a linear-with-time increase of the mean 
kinetic energy. In this case there is very little to distinguish 
resonant from antiresonant evolution lH |9l [lOl [TT] . Narrower 
initial momentum distributions are achievable, however For 
example: the Rubidium 87 Bose-Einstein condensed sample 
of Duff^y et al |27| had an initial width of w ~ 0.043ft/: 
(with A - An/K - 780nm; ID Raman-cooled caesium with 
w ~ 0.062ft/r ll62l . and Bose-Einstein condensed caesium 
with ~ 0.008ft/: (W\(A = 895nm) is achievable; and sam- 
ples of ID Raman-cooled sodium with w ~ 0.Q0l6fiK 1 63] 
(A = 589nm) have been reported. Experimental exploration 
of the regimes described in this paper is therefore quite possi- 
ble. 
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The simplest approach is to use that 

F(J3)\k J dxF(^)\x){x\k + /?>. 



(Al) 



Substituting the Floquet operator (20 1 into Eq. (All yields 

F(fi)\k + P) ^ e-'"'^''^ J dxe-''''^^e'^-''°''^'^^\x){x\k+li), (A2) 



where we have defined 



Ky = 7i(\ +2/3)e. 



(A3) 



The kicking term in the Floquet operator introduces a spatially 
dependent phase. 



e 



i4>d cos(Kx) 



x) — \x}e' 



i^ti cos{/lx) 



(A4) 



and the free evolution leads to a position kick ||74|| . i.e., 

e-''^'^y\x) = e'^'^''|x + y). (A5) 



The matrix element {x\k + /3) in Eq. ( A5 i is known [see 



Eq. (26 1], so we find that 




27T J 

(A6) 



2. Application of subsequent Floquet operators 



The time evolution of state (jASj due to a second d-kick, and 
free evolution for time T - {Tt/2, can be found in a similar 
way, to be 




F^mk + yS> = y --e'2/5^r^-«./^^^ J ^_,|^ + 2y) ^^^^ 

^ ^i(k+l3)K.x^iM<:os{Kx)+cos(Kx+Ky)] 

Subsequently, the effect of the full time-evolution operator 
F"i/3) is 



F"mk+/3}^^l^e'"'"'ye 



'""P"' ^ dx\x + nyye'^'^P'"''' 



(A8) 



APPENDIX A: EVALUATION OF THE FLOQUET 
OPERATOR MATRIX ELEMENT 



1. Application of one Floquet operator 

In this Appendix we calculate the matrix element 
{x\F"(J3)\k+/3}. We begin by considering the effect of applying 
the Floquet operator F(J5) to the momentum eigenket \k + j3). 



3. Spatial representation 



The final state (A8i can be projected onto the ket \x'), and 



the orthogonality relation {x\x') - 6{x - x') means that the 
integral can be evaluated. That yields 



K ,,, 



{x\F"(B)\k + B}^J ^e'>'/^^frg-""/i^fg'(«:+/3)^:(A-«y) 

27T (A9) 



X e 



itf>j^cos(Kx-Kny)+X"jJi cos(Kx-nK-y+ jK-y)^ 
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Equation ( A9 1 can be simplified by defining the variable 

q. = (n- j)Ky, (AlO) 

to give 



MI3Ky -irmf t i(k+l3)(Kx-qo) 



V 2n 



(All) 



Defining f = 2"=d cos qj and ^ = 2"=o sin^j, Eq. (All 
be further simplified to 



{x\F"(J3)\k + y8) = y ^e'"'^Ky^-mn,3 t ^i{k^l3)(Kx-q„) ^^^^^ 
^ ^i(l>j^cos(Kx)+i(l>j(sm{Kx) 



4. Bessel function expansion 

Invoking the Bessel function expansions 



]=- 



(A13) 



j=- 



we can rewrite Eq. ( A12 1 as 



j=- 



2 



Using Graf's addition theorem ll64l . Eq. ( |A14| l becomes 



{x\F"mk + B)^J —e'"'^'^ye-"""^ e ^i(k+p)(Kx-q„) 



(A15) 



where oj andx are real and defined by we'-^ = (pd(i^ + 0- It is 
useful to define the function 



n-l 
i=0 



(A16) 



from which ^ and ^ can be deduced via ^ = (ji+fi*)/2 and f : 



(// - A'*)/2/. Substituting [from Eq. (AlOi] into Eq. (A16i 
we find that 



(A17) 



and evaluating the geometric sum gives 

1 - e-'"'^^ 



inKy 



1 - e-'^r 

Using standard trigonometric identities, we find that 

^ = sin(nT)[cos(nT) cot T - sin(nT)], 
^ = sin(«T)[cos(«Y) + cotTsin(«Y)], 



(A18) 



(A19) 



where we have defined the parameter T = Kyl2 [this is con- 
sistent with the definition of Eq. ([28]l]. Finally, 



<l>d 



sinT 



(A20) 



and substituing back into Eq. (A15l, the matrix element 
{x\F"(J3)\k + y8) is determined to be 



K 



2n 



ilnpr -hmfj-l i(k+l3)(Kx-2nr) 



Z°° ::fc^ I sin(nY)\ ,. 



-'■>(«+ 1)T 

(A21) 



APPENDIX B: MOMENTUM-MOMENT EVOLUTION FOR 
A LOW TEMPERATURE ATOM CLOUD 

1. Second-order momentum moment evolution for a 
momentum eigenstate with k = Q 

We consider a system that is initially prepared in the mo- 
mentum eigenstate ^), i.e., k - Q and y6 e [-1/2, 1/2). After 
a time riT, the second-order moment is given by [see Eq. pT| ) 
with q - 2 and ^ = 0] 



(A14) {p\^n^K^Y. 



Jj\(pd 



sin(nY) 
sinY 



(Bl) 



where Y is defined in Eq. (28 i. For compactness, we define 
the p dependent variable 



sin(nY) 
sinY ' 



(B2) 



and then Eq. (Bl i can be expanded to give 



<p2)„ = n^K^ £ J](j^)[f- + 2jl3+li\ (B3) 



j=- 



To evaluate Eq. ( B3 i, we consider each the three terms on 
the right-hand side, in reverse order To evaluate the third term 
we use the identities [64J 



and 



7_/77) = i-iyjjiJ]), 

oo 

jI{ii) + 2Yj]{iT)^1, 



(B4) 



(B5) 
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to find that 



(B6) 



The second term can be evaluated, using Eq. (B4i, to give 

CO 

^ Jj(r])2jfi = 0. (B7) 

To evaluate the first term in Eq. (B3 1 we use the identity Ii64l 

(B8) 

1 

and Eq. (|B6|l, to find that 



(B9) 



2j 

Jj+iiif) + Jj-i{ri) = —Jjiji), 



=— oo 

Using the addition theorem ll64l 

oo 



(BIO) 



along with the identity Jji-rf) = (-iyjj(rj), and Eq. (B4i, 
Eq. (B9 1 becomes 



X fjjiri) = [1 + ^2(0)] = -/7^ (BID 



Therefore, the second-order momentum moment evolves as 



f>', sin2(„T) 



2 sin^Y 



(B12) 



2. Second-order momentum moment evolution for a square 
distribution 



be evaluated as shown in Sec. B 1 to yield 



(e^ r 



4e sin^T , 



(B15) 



To evaluate the remaining integral, we make a change of vari- 
ables such that 



A. f 



..«l+2.)/2 ^j„2(„Y) ^ 



3 Aent J„n^_2e)ii sin^T 
Then using the expansion 1551 

,2/ >v.s 1-1 



sin''(nY) 
sin^T 



= n -H 2 ^(n — ^m) cos 2mT, 



(B16) 



(B17) 



where m is integer, we find that 



1 



i2 n-1 



2e7r/' 



(-1)"" sm(2mn{e) 

m I 



(B18) 



3. Second-order momentum moment evolution for a Gaussian 
distribution 

Here we consider the evolution of the second-order momen- 
tum moment for the case where the initial momentum distri- 
bution is Gaussian (see Sec. |llB4| . In particular, we consider 
an ultracold gas for which the momentum width of the initial 
distribution is very small compared to the two-photon recoil 
momentum tiK, i.e., w <K 1. Thus, initially the majority of 
atoms have k = 0, and their /3 values are typically much less 
than the maximum possible value \/3\ -1/2. The second-order 



momentum moment evolution is described by Eq. ( 37 i with 



'V2^ 



(T])e 



-{j+py. (B19) 



We consider a square initial atomic momentum distrubition 
of the form 



0, 



> e. 



(B13) 



The initial standard deviation is e/ V3, and the second-order 
momentum moment evolves as 



2 1-2 



Xe uo 
dl5Y,j]^r^){j+fif, (B14) 



2e 



where we have taken Eq. ( 35 1 with q - 2 and ^ = 0, and with 



DkiP) given by Eq. (B13 i. The sum over j in Eq. (B14i can 



In the regime where w <k 1, we can approximate Eq. (B19l 
by extending the limits of the p integral to infinity, since the 
additional contribution to the integral will be negligible. Eval- 



uating the sum over j as shown in Sec. B 1 yields 



w 



V2^ 



'2 



2i-2j,2 



Zw ' — ( 



i r,^w),-.v2.. 

n 



(B20) 



sin^T 



In the case that an antiresonance occurs at y6 = 0, i.e., £ is 
odd, we can write T = 7r(m +1/2-1- /3£), where m is an integer 



[see Eq. (|4T])]. Thus, sinY = (-1)'" cos(7rjeO and sin(«Y) = 
(-1)"'" sin(«7r/2 + hk/S^). In the regime where <k 1 /2, we 
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can write that sin T « (-l)™, which leads to the considerable 
simplification that 



w 



Zw 



oo 



2 , 'fl 



1 



1 p 

V27r J-cx 



t/yS cos(2nT)e 



(B21) 



It is reasonable to take \f^{\ <s 1 /2 since w <k 1 implies that the 
only nonzero contributions to the momentum moment evolu- 
tion are from small 

Using that T - 7i(m +1/2 + I3£), we write cos(2»T) - 



cos(«7r) cos{2njT/3{), and then the integral in Eq. (B21 1 can be 
evaluated as follows: 



) a; W H 



— w + 



— w + 



1 - 



cos(n7r) 
w V2^ 
cos(n7r) 



w 



V2^ 



CO 

Xoo 
OO 



(B22) 
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